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Abstract: We present a study of showers initiated by electrons, pions, kaons, and protons with 
momenta from 15 GeV to 150 GeV in the highly granular CALICE scintillator-tungsten analogue 
hadronic calorimeter. The data were recorded at the CERN Super Proton Synchrotron in 2011. 
The analysis includes measurements of the calorimeter response to each particle type as well as 
measurements of the energy resolution and studies of the longitudinal and radial shower develop¬ 
ment for selected particles. The results are compared to Geant4 simulations (version 9.6.p02). In 
the study of the energy resolution we include previously published data with beam momenta from 
1 GeV to 10 GeV recorded at the CERN Proton Synchr'otron in 2010. 

Keywords: Calorimeter methods; Detector modelling and simulations I (interaction of radiation 
with matter, interaction of photons with matter, interaction of hadrons with matter, etc); Particle 
identification methods. 
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1. Introduction 


A sampling hadronic calorimeter with highly granular readout, as required for Particle Flow Anal¬ 
ysis (PFA) OJ, 0, and with tungsten as absorber, is considered for experiments at future multi-TeV 
e+e^ colliders such as the Compact Linear Collider (CLIC) |0]- In these experiments the use of 
tungsten as dense absorber provides a calorimeter system sufficiently deep to contain jets of TeV 
energy, with a diameter fitting inside a solenoid of acceptable size. 

This paper presents the analysis of data obtained in test beam experiments of the CALICE tung¬ 
sten analogue hadronic calorimeter prototype (W-AHCAL) 0, EL 01 using mixed beams containing 
electrons, muons, pions, kaons, and protons with a momentum range from 15 GeV to 150 GeV (in 
this paper, the natural system of units with /i = c = 1 is used). The analysis includes measurements 
of the calorimeter response to each particle type, the energy resolution, and studies of the longi¬ 
tudinal and radial shower development. A purpose of the analysis is the comparison of data with 
Geant4 HZlil simulation models. Based on the comparisons, the reliability of the models used in 
detector studies for future collider experiments can be assessed and the simulation models can be 
further improved. 

Section ^ describes the experimental setup used in the test beam experiment. The calibration proce¬ 
dure is discussed in section 0 and the data selection in section ^ The Monte Carlo (MC) simulation 
used in the comparison with data is introduced in section ^ Section ^ lists the systematic uncertain¬ 
ties relevant for the studied observables. The analysis results for positrons are shown in section |7| 
and for hadrons in section ^ In section we compare the W-AHCAL response for different 
particle types, and in section we draw conclusions. 

2. Experimental setup 

The W-AHCAL setup comprises a stack of 38 layers of absorber plates interleaved with 0.5 cm 
thick scintillator tiles, read out by wavelength shifting fibres coupled fo silicon phofo-mulfipliers 
(SiPMs). Each absorber plate is 1 cm fhick and is made of a fungsfen alloy consisfing of 92.99% 
fungslen, 5.25% nickel, and 1.76% copper, wifh a densify of 17.8 g/cm^. The nuclear inferacfion 
lengfh of fhis alloy is Ai = 10.80cm and fhe radiafion lengfh is Xq = 0.39cm. The scinfillafor 
files of each layer are mounted info sfeel casseffes, wifh 0.2 cm fhick walls and are connecfed fo 
prinfed circuif boards (PCB). The overall fhickness of a calorimefer layer is approximately 2.5 cm. 
The combined sfrucfure of absorber, sfeel supporf, PCB, and scinfillafors has an effective infer- 
acfion lengfh of Aj = 19.0 cm and radiation lengfh of Xq = 0.88 cm such fhaf one calorimefer 
layer corresponds in fofal fo 0.13 X\ and fo 2.8 Xq. The overall dimensions of fhe profofype are 
0.9 X 0.9 X 0.95 m^, amounting fo 4.9 Ai and fo 108 Aq. The high granularify of fhe detector is 
ensured by the 3x3 cm^ tiles placed in the centre of each active plane, surrounded by 6 x 6cm^ 
and 12 x 12cm^ tiles at the edges. The last 8 layers of the stack are equipped with 6 x 6cm^ tiles 
in the centre surrounded by 12 x 12 cm^ tiles. In total the calorimeter is equipped with 7608 cells. 
The temperature of each layer is monitored by five sensors in order to correct for the temperature 
dependence of the SiPM response. The response of all SiPMs is calibrated and monitored with a 
EED system 0]. 

The data were taken in 2011 at the H8 beam line |@, HO] of the CERN Super Proton Synchrotron 
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Figure 1: Sketch of the CERN 2011 test beam line (not to scale), where Sc stands for scintillator 
and WChfor wire chamber. The beam enters from the right. 


(SPS). The primary proton beam from the SPS hits a target 605 m upstream of the W-AHCAL 
prototype. The secondary beam can then either be used for the experiment or it can hit a secondary 
target 477 m upstream of the W-AHCAL prototype. A momentum selection and focusing system is 
used to deliver the secondary or tertiary beam to the experiment. Either beams of e* with momenta 
between 10 and 40 GeV or mixed beams of tt*, K^, and protons between 10 and 300 GeV were 
used. 

The data analysis is restricted to e+ events with beam momenta from 15 GeV to 40 GeV as well as 
to positively-charged hadron events with beam momenta from 25 GeV to 150 GeV for which high 
selection purity can be achieved and leakage effects are limited (cf. section Negative-polarity 
data have smaller statistics due to less allocated beam time, but they are included if necessary for 
the study of detector-related systematics. 

A sketch of the CERN SPS test beam setup is presented in figure |T]. A right handed coordinate 
system is used, with the z-axis given by the beam axis, the x-axis pointing horizontally and the 
y-axis pointing up. The beam passes through two Cherenkov threshold counters, two trigger scin¬ 
tillators and a tracking system of three delay wire chambers flTTlI before reaching the W-AHCAL. 
The Cherenkov counters are filled with He gas. The Cherenkov radiation is detected by photo¬ 
multiplier tubes that are operated with discriminators at fixed thresholds. The Cherenkov signals 
are used offline for hadron idenfification, one Cherenkov counter being used for the separation of 
pions from kaons and the other one for the separation of kaon from protons. The beam trigger is 
defined by the coincidence of two 10 x 10 x 1 cm^ scintillator counters. The information from the 
three 11x11 cm^ wire chambers is used offline to reconstruct the trajectory of the incident particle 
and predict its position on the calorimeter surface. 

In the analysis of the W-AHCAL energy resolution, data from the 2010 measurement campaign Unil 
at the CERN Proton Synchrotron (PS) with beam momenta from 1 GeV to 10 GeV are also in¬ 
cluded. The experimental setup for these low energy measurements was slightly different. The 
calorimeter was equipped only with the first 30 layers accounting for 3.9 Aj and 85 Aq. The 
Cherenkov counters were placed at 3.5 m and 7 m in front of the experimental setup and filled 
with CO 2 to better distinguish the lower momentum particles. 
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3. Calibration 


The response of all calorimeter cells is equalised to a common physics signal produced by mini¬ 
mum ionising particles (MIP), which is used throughout the paper as the energy scale, being the 
most natural one for the comparison with Monte Carlo simulations. For this calibration, data from 
dedicated muon runs are used, as muons are considered to be minimum ionising particles. Several 
steps are necessary to translate signals recorded by the SiPM measured in ADC counts into de¬ 
posited energy in units of MIP. 

First a track finding algorithm [ffUl is used to reconstruct the muon track. The energy distribution 
in each cell for hits associated to a reconstructed muon track is fitted with a Landau function con¬ 
volved with a Gaussian function. The MIP calibration factor of a single cell i, [ADC] is given 
by the most probable value of the signal of a muon in that cell, expressed in ADC counts. The 
calibration of a single cell i is then performed according to 

•/-p(A,[pixel]), (3.1) 

where 


• A, [ADC] is the pedestal subtracted amplitude registered in cell i, in units of ADC counts, 

• A^^^[ADC] is the MIP amplitude in cell i in ADC counts, called the MIP calibration factor 
and extracted from the dedicated muon runs, 

• /resp(A,[pixel]) is the SiPM response function which corrects for the non-linearity of the 
SiPM response due to the limited number of pixels (1156) and the finite pixel recovery time 
(20 — 500 ns). This function acts on the amplitude expressed in number of pixels, and returns 
the saturation correction factor to linearise the amplitude in MIPs. Its asymptotic behaviour 
is described using a double-exponential fit function and a linearised extrapolation for large 
values nnil . 


• The amplitude in pixels is obtained by dividing the amplitude of a cell by the corresponding 
SiPM gain Gj [ADC/pixel] 


A,-[pixel] 


A,-[ADC] 

G; [ADC/pixel]' 


(3.2) 


The gain values G,- [ADC/pixel] are obtained from fits of photo-electron spectra taken with 
low intensity LED light provided by the calibration and monitoring LED system HnH . 


In the following data analysis, only cells with energy depositions above a threshold of 0.5 MIP are 
considered. 


3.1 Saturation scaling factor of the SiPM response function 

The response of all SiPMs was measured before mounting the sensors on the scintillator tiles [0]. 
Due to geometric alignment effects, the effective maximum number of pixels for mounted SiPMs 
is on average 20% lower than the value for bare SiPMs HnH . Eor this reason, the response functions 
which correct for saturation in the SiPM response of all AHCAE cells are scaled by a factor s which 
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Energy sum (MIP) 


Figure 2: Distribution of the energy sum over 
all cells for 40 GeV data. Data reconstructed 
using individual and average saturation scaling 
factors are compared to simulations. 



Figure 3: Hit energy distribution for the most 
central cell with respect to the beam position in 
layer 3 for 40 GeV Data reconstructed using 
individual and average saturation scaling fac¬ 
tors are compared to simulations. 


has been determined in situ using the LED calibration system. While for steel and low energy 
tungsten data the usage of an average factor of 5 = 0.8 was found to be sufficient, the high energy 
tungsten data exhibit higher single cell energies, such that the usage of the exact scaling values is 
required for a correct reconstruction. Hence, individual scale factors are used where available (in 
60% of the cells) and the average value of 5 = 0.8 otherwise. The impact of using individual or 
average scale factors is largest for the highest positron energy, as shown for the total energy sum 
distribution of the 40 GeV e+ data in figure The dafa evenf selecfion and fhe simulafion used in 
fhis sfudy are infroduced in secfions ^ and The effecf is particularly sfriking in single-cell hif 
energy disfribufions of cells which have an individual scaling factor fhaf differs significanlly from 
fhe average scaling faclor. As an example, fhe hif energy disfribufion for 40 GeV e+ for fhe mosf 
cenfral cell wifh respecf to fhe beam posifion in layer 3, as evaluated using fhe average (0.8) and 
individual (0.91) scaling factor is shown in figure ^ 

Differences befween fhe shower properties reconsfrucfed using individual and average safurafion 
scaling faclors decreases fowards lower beam momenfa and are smaller for hadron showers fhan 
for elecfromagnefic showers. For fhis reason, in fhe reconsfrucfion of fhe PS dafa sefs wifh beam 
momenfa befween 1 GeV and 10 GeV used in fhe sfudy of fhe energy resolution in secfions 0 and 
1^ average scaling values as used in IIT^ were adopfed. In fhis momenfum range, fhe visible energy 
and fhe energy resolution obfained using individual and average scaling agree wifh each ofher on 
fhe sub-percenf level. 

3.2 Temperature correction 

The temperature inside the calorimeter is monitored using five sensors per layer. Details about 
the temperature measurements are given in [1T3I] . Temperature variations of several degrees were 
observed with time and along the beam direction. Due to the temperature dependence of the SiPM 
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response, the MIP calibration factors need to be corrected for temperature variations. The temper¬ 
ature dependence of the MIP calibration factor of a given cell i can be described by the formula 


AfP(rrun) = AfP(ro) +C- (Trun - To), 


(3.3) 


where 


• A“^^(rnin) is the MIP amplitude in cell i in ADC counts at the temperature Trun, 

• A“^^(ro) is the MIP amplitude in cell i in ADC counts, obtained at the reference temperature 
To, i.e. at the average temperature during the muon runs of the sensor closest to the cell, 

• Trun is the temperature of the sensor closest to the given cell in the given run in °C, 

• C is the MIP temperature coefficient determined for the layer of the given cell in ADC/°C. 

Due to limited statistics, the MIP temperature coefficients are not determined for each cell indepen¬ 
dently but estimated in average for all cells in one calorimeter layer, using the method described 
in m. For a few calorimeter layers, no reliable estimations for the MIP temperature coefficient 
could be determined. For these specific layers, the average of all other coefficients is used. The 
average MIP temperature coefficient considering all layers is —4.5%/°C before the temperature 
correction, and better than 0.06%/°C after. A similar procedure is used to correct the temperature 
dependence of the gain calibration coefficients. 

4. Data selection 

The data analysis was performed for e+ events with beam momenta from 15 GeV to 40 GeV, and 
proton events with beam momenta from 25 GeV to 150 GeV, and K+ events at 50 GeV and 60 GeV. 
For these data sets, high selection purity can be achieved [HBl and leakage effects are limited (cf. 
section 

For most of the particles and beam momenta, 20000-100000 events (figure ^ pass the selection 
steps introduced in the following. For 60 GeV and 80 GeV, high statistics hadron data sets were 
recorded with approximately up to 700 000 events after event selection. Beam momenta with par¬ 
ticularly low number of events after event selection are 40 GeV e+ (6000 events), 50 GeV K+ 
(3 000 events), and 150 GeV (5 000 events). 

4.1 Electromagnetic showers 

Positron events were recorded at the SPS with a dedicated positron beam that has a small contam¬ 
ination from muons. To increase the purity of the e+ data, an optimised cluster reconstruction and 
subsequent event selection procedure is pursued. 

Three-dimensional clusters are identified using the nearest-neighbour algorithm developed in flTTIl . 
Cells with an energy above a threshold of 2 MIP are used as seeds for the clustering. All active 
cells which surround these seed cells and have an energy deposit above a threshold of 0.5 MIP are 
added to the cluster. 

The reconstructed e+ event candidates are selected if they fulfil the following requirements. First, 
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Figure 4: Number of events after selection for a given particle type. 


the number of identified clusters in the whole calorimeter is required to be one. Electromagnetic 
showers form a cluster near the calorimeter front face, hence in the second step events are rejected 
if their centre-of-gravity in the z-direction (Zcog) exceeds 160 mm, which rejects muon events. Here, 
Zcog is defined as 


Zcog — 


L Mlits P ’ 
/=! 


(4.1) 


where Ei is fhe visible energy deposifed in fhe i-th acfive calorimefer cell, Zi is ifs position in fhe 
reference frame of fhe calorimefer wifh fhe W-AHCAL fronf plane af z = 0 mm, and AhUs is the 
number of W-AHCAL cells which have an energy deposition above the threshold of 0.5 MIR In a 
third selection step, the events are required not to contain an identified frack which could indicafe 
a muon or a hadron evenf in which eifher fhe muon or fhe hadron fhemselves or particles inside 
fhe hadron shower fraverse (parfs of) fhe calorimefer sfack as a minimum ionising parficle. As a 
lasf requiremenf, Ahits should be inside fhe range defined in fable |T]. This is done in order fo rejecf 
fhe few evenfs in which eifher fhe election shower is inifiafed before enfering fhe W-AHCAL or 
evenfs fhaf are particularly noisy. A fypical e+ evenf affer all selection cufs is shown in figure || 
Depending on fhe beam momenfum, befween 60% and 80% of all position candidafe evenfs pass 
fhe described position selecfion. 

To reduce fhe influence of noise in fhe selecfed e+ evenfs, only calorimeter cells wifhin fhe firsf 20 
calorimefer layers (corresponding fo 56 Xq) and wifhin fhe cenfral 10 x 10 files of 3 x 3 cm^ size 
(corresponding fo more fhan 3.5 effecfive Moliere radii) are used fo calculafe fhe position shower 
properties discussed below. This sfrafegy is appropriafe because e+ showers in fhe energy range 
sfudied here deposif energy mosfly wifhin fhe firsf 5 fo 10 calorimefer layers, fhe beam is cenfred 
on fhe calorimefer cenfre, and fhe widfh of fhe beam profile is nof more fhan fhree files. 
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Table 1: Cuts on the minimum and the maximum number of hits applied for the selection. 
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Figure 5: Display of a 15 GeV e^ event. 
The colours of the W-AHCAL cells indicate 
the deposited energy, red being the highest 
(more than 5.4 MIP). 



Figure 6: Display of a 80 GeV 7t+ event. 
The colours of the W-AHCAL cells indicate 
the deposited energy, red being the highest 
(more than 5.4 MIP). 


4.2 Hadronic showers 

Hadron events are recorded at the SPS with a mixed beam containing mostly K+, protons, 
and |r+ from in-flight decays. Contamination from e+ was effectively suppressed by using a lead 
absorber in the beam line. The absorber had a thickness of 8 mm for the 25 GeV and 30 GeV data 
sets and 18 mm for higher beam momenta. Even for an absorber thickness of 8 mm the effect of 
the resulting e+ impurity on the shower observables is negligible. Contamination from p+ from 
in-flight decays is rejected based on information from the W-AHCAL itself, as described below. 
The first step in the hadron selection is based on particle identification information obtained from 
the two Cherenkov threshold counters installed upstream of the W-AHCAL prototype as described 
in section ^ A detailed description of the Cherenkov selection strategy and results can be found 
elsewhere ma. The event samples of pions, kaons, and protons selected based on the Cherenkov 
threshold counters contain between 85% and 100% of the selected hadron type. 

Events are then rejected if they do not contain a reconstructed cluster (cf. section O)- After that, 
the layer of primary interaction of the reconstructed shower, i.e. the shower start, introduced in 
UTEH . must be in the first three layers of the calorimeter, except when measuring the longitudinal 

















profiles where an alternative selection is used (see section |8.1.2| ). The latter requirement allows 
for a selection of hadron showers that are mostly contained in the W-AHCAL stack as discussed 
below. Both requirements reject muon events. The remaining muon contamination is reduced by 
exploiting the fact that muons show a different signature in the highly granular W-AHCAL from 
hadron events, i.e. they deposit energy only in comparatively few cells of the full W-AHCAL stack 
and the mean depth of all energy depositions given by the centre-of-gravity of all hits in the z- 
direction has a different distribution from that for hadrons. In the (Ahits : Zcog)-plane, all events 
below a threshold given by Ahits + ci{pheam) *Zcog < (^(Pbeam) (where a and b depend linearly on 
Phsam) are rejected. While removing few remaining muon events, the cut does not have any impact 
on the central part of the energy sum distribution of the hadron data. An example 7C+ event after all 
selection cuts is shown in figure Depending on the beam momentum, between 5% and 30% of 
all hadron candidate events pass the described hadron selection, mostly influenced by the variation 
of the beam’s muon content. 

The W-AHCAL prototypes used in the two test beam campaigns at the PS and SPS, with respective 
total lengths of 3.9 Aj and 4.9 Ai, do not fully contain all hadronic showers at the studied beam 
momentum ranges of 1-10 GeV and 15-150 GeV, even if the showers start within the first three 
W-AHCAL layers. The remaining longitudinal leakage of energy after applying the shower start 
cut has been studied using simulated hadron showers in a long W-AHCAL prototype of 100 layers 
corresponding to 13 Aj. The energy fraction contained in 30 or 38 layers of this long prototype is 
shown in figure |7| for simulated proton showers. Events are selected by requiring that the shower 
starts within the first three layers of the calorimeter, as also required in the analysis of hadron data. 
Similar results are found for different Geant4 physics lists and hadron types. 

As can be seen in figure 0 for protons, 5.5% of the shower energy measured in 100 layers is expected 
to leak out of a 38 layers calorimeter at 150 GeV. In order not to be dominated by leakage effects in 
the characterisation of the W-AHCAL, we limit the momentum range studied to this upper value of 
150 GeV. We estimated that leakage has a L0%-5.5% effect on the W-AHCAL response (figure [^l, 
a 0%-17% effect on the energy resolution for the studied momentum range and it has a negligible 
effect on shower shapes. The detector simulations introduced in section 0 use the same W-AHCAL 
geometry as in data such that leakage effects are taken into account. 

4.3 Combination of data recorded at the CERN PS and SPS 

To better constrain the W-AHCAL energy resolution fit results discussed in sections 0 and the 
W-AHCAL data of the test beam campaign at the CERN PS with beam momenta from 1 GeV to 
10 GeV nnH are included. Eor this, the PS data were selected based on the SPS selection cuts 
described above. The PS analysis results obtained in this way are slightly different from but con¬ 
sistent within the uncertainties to the previously published results flO . In the combination of the 
PS and SPS data for the energy resolution, the difference in the energy resolution between 30 and 
38 layers has been taken into account by introducing a systematic uncertainty of 2% on the PS data 
points. This uncertainty was estimated using hadron shower simulations at PS energies as intro¬ 
duced in section by comparing the energy resolution of a W-AHCAE of 30 and 38 layers. Eor 
SPS energies we observe that the difference in the energy resolution obtained with 30 layers and 
38 layers in data is well described by MC. 


-9- 




Geant4 simulation 



Figure 7: Energy fraction detected versus beam momentum of simulated proton (p) showers within 
30 (38) layers of the 100-layer deep W-AHCAL at PS (SPS) energies and energy fraction detected 
within 30 layers of 38-layer deep W-AHCAL at PS energies. The showers are required to start 
within the first three layers of the calorimeter. 

5. Monte Carlo simulation 

The simulation includes the W-AHCAL, the two trigger scintillators and the three wire chambers. 
Particles of a given type and energy are generated at z = —55 m. The transverse beam profiles in 
the simulations are chosen such that they reflect the beam profiles measured in the corresponding 
data runs. Interactions are modelled using the QGSP_BERT_HP and the FTFP_BERT_HP physics 
lists of Geant4 version 9.6.p02. These lists are described in llOl . The extension HP indicates that 
the physics lists use a data driven high precision neutron package for the transport of neutrons be¬ 
low 20 MeV down to thermal energies. This is crucial for the simulation of the high-A absorber 
tungsten where many more neutrons due to spallation are produced than in steel absorbers llIPIl . 
The size of the Monte Carlo data sets are chosen such that they reflect the approximate number of 
events of the individual data sets used in the analysis. The same event selection as introduced in 
section ^ for data is applied to the simulated data sets; similar selection efficiencies in simulation 
and data are observed. 

For a valid comparison between simulation and data, realistic detector effects need to be consid¬ 
ered at the generation and the digitisation level. In the generation step, the signal shaping time of 
the readout electronics is emulated in the Monte Carlo simulation by accepting only hits generated 
within a time window of 150 ns llTTII . In addition, saturation effects in plastic scintillators, described 
by Birks’ law, are applied lEOl . In the digitisation, the same detector granularity, calibration values 
and dead or uncalibrated cells are considered as for the reconstruction of the experimental data. 
The simulated energy is converted into MIP based on a conversion factor estimated in muon simu¬ 
lations imH . Fight leakage between neighbouring scintillator tiles is also taken into account [ETll . 
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During data taking, some events were recorded with random triggers in order to monitor the noise 
level in the calorimeter HHD- This noise is not subtracted from the data but it is accounted for by 
superimposing random trigger events from the relevant run period on simulations. 

In contrast to previous AHCAL studies, in this study an additional treatment of saturation effects 
is included in the simulations. In the data reconstruction, the asymptotic exponential form of the 
SiPM response function introduced in section || is extrapolated with a linear function for very large 
signal amplitudes in order to avoid unphysical results. For the large amplitudes occurring when 
using tungsten absorbers and at large beam momenta, this potentially leads to an under-correction 
of saturation effects in data. In previous AHCAL studies, the same form of the SiPM response 
function was used in simulation and reconstruction such that the effect cancels out in the case of 
MC, but not for data. By retaining the asymptotic form in the simulation while using the linearised 
function in the reconstruction, saturation effects observed in data are better reproduced by the sim¬ 
ulations. 

6. Systematic uncertainties 

The most relevant sources of systematic uncertainties for the W-AHCAL test beam data are the 
limited knowledge of the saturation scaling factor of the SiPM response function, the uncertainty in 
the MIP calibration, the stability of the W-AHCAL response over time, and the choice of the event 
selection based on the shower start layer. For the Geant4 simulations, the most relevant systematic 
uncertainty originates from the simulation of the cross-talk between neighbouring scintillator tiles. 
The systematic uncertainties due to all relevant sources and all discussed observables are studied 
and quantified. In the following, for brevity, only the systematic uncertainty on the reconstructed 
visible energy and, in specific cases, on the longitudinal or the radial energy profile are quoted. 

6.1 Saturation scaling factor of the SiPM response function 

For electromagnetic showers in the W-AHCAL, due to the dense absorber with about 3 radiation 
lengths per layer, most of the shower energy is deposited in a few cells in the first calorimeter lay¬ 
ers. Hence, uncertainties in the SiPM response function are expected to have significant impact on 
the reconstructed e+ shower properties. 

As discussed in section ^7T| , the response functions for all SiPMs were measured before mount¬ 
ing the sensors on the scintillator tiles. These are then scaled with a saturation scaling factor s, 
accounting for the imperfect mounting of the SiPMs to the wavelength shifting fibres of the scin¬ 
tillator tiles. The distribution of the measured saturation scaling factors with a mean of 5 = 0.8 has 
a standard deviation of 0.09, reflecting cell-to-cell variations and measurement uncertainties. Indi¬ 
vidually measured values are available for 60% of the cells. The uncertainties on these factors are 
taken to be ±0.5 standard deviations of their distribution (±0.045) for the individually measured 
factors and ±1 standard deviation (±0.09) for the saturation scaling factors for the unmeasured 
cells. 

Systematic uncertainties in the reconstructed deposited energy are estimated by reconstructing the 
data using scaling factors varied by the full range of their respective uncertainties, denoted as high 
and low scaling. An increase in the scaling factor results in a decrease of the reconstructed energy 
and vice versa. 
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Figure 8: Comparison of the longitudinal energy profile for positron-induced showers at 40GeV 
(left) and pion-induced showers at 100 GeV (right) when using SiPM scaling factors varied within 
the full range of their uncertainties. The bottom panels (same for figures ^ andfT^fT^) show 
the ratio between the distributions for the different choices of scaling factor. 


Figure shows how the maximal modification of the saturation scaling factors within their uncer¬ 
tainties increases or decreases the reconstructed energy deposition per layer in longitudinal energy 
profiles for e+ (left) and 7t+ (right) induced showers. The observed variations, S^vis per iayer> are 
treated as the limits of a rectangular distribution and the resulting uncertainties are used as sys¬ 
tematic uncertainties on the longitudinal energy profile. As an example, the observed variations 
and the resulting systematic uncertainty on the layer-by-layer energy deposition for e+ showers for 
the layers with the maximal variation are listed in table The systematic uncertainty increases 
as a function of the beam momentum, since saturation of single cells is reached sooner at higher 
momenta and saturation scaling factors play an increasingly important role in the cell-energy recon¬ 
struction. The layer-by-layer uncertainties in hadron showers at the same beam momentum have 
lower values than in the denser e+ showers as more cells contribute to the energy sum of one layer 
and saturation is less common. For hadrons, the maximal layer-by-layer systematic uncertainties 
range from 1.0% at 25 GeV to 3.5% at 150GeV. 

In the study of the detector response discussed in sections 0 and 0, the average visible energy in a 
particle shower measured with the W-AHCAL, (£'vis)^ for particle showers of one particle type and 
beam momentum is obtained by fitting the corresponding energy sum distribution to a Gaussian 
function. The event’s energy sum is obtained by summing up all cell energies for all cells i 
which have an energy above threshold. 

The systematic uncertainty on the average visible energy is approximated by the quadratic 
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Table 2: Maximal impact of the SiPM saturation scaling factor modification on the longitudinal 
energy profiles for e^ showers and the resulting maximal systematic uncertainty on the layer-by- 
layer energy deposition. 


Pbeam (GcV) 

^^vis per layer, max (%) 

Syst. uncertainty (%) 

Low scaling 

High scaling 


15 

+3.6 

-2.8 

1.8 

20 

+4.8 

-3.7 

2.5 

30 

+6.6 

-5.4 

3.5 

40 

+7.5 

-6.1 

3.9 


Table 3: Systematic uncertainty on the average visible energy of the full showers due to the 
uncertainty on the SiPM saturation scaling factor. 


Pbeam (GcV) 

se 

(%) 

Low scaling 

High scaling 

15 

+ 1.4 

-1.1 

20 

+ 1.9 

-1.4 

30 

+2.7 

-2.1 

40 

+3.0 

-2.5 


sum of the layer-by-layer response variation given by 




default ^£,modified)^ 


(6.1) 


where ^ is the index of the layer, is the visible energy in layer i as reconstructed with default 
or modified saturation scaling values in all cells and ^^ax is 20 for e+ and 38 for hadrons. The 
resulting systematic uncertainties on the visible energy for the analysed e+ energies are given in 
table ^ For the visible energy of hadron showers, in which many more cells contribute to the full 
shower energy, the impact is reduced and ranges from 0.4% at 25 GeV up to 1.5% for 150 GeV. 


6.2 MIP calibration 

The MIP calibration factors have two types of uncertainties: a statistical uncertainty on the mea¬ 
surement itself and a systematic uncertainty associated with the method. The statistical uncertain¬ 
ties are uncorrelated between the different cells. In 111411 a MIP systematic uncertainty of 2% was 
found, due to imperfections of the parametrisation of the MIP line shape and other effects such as 
binning or bias in the muon track selection. This uncertainty affects all cells in the same way. 
Muon calibration runs, which are used for a cell-wise MIP calibration as introduced in section 0, 
were recorded twice, in summer and autumn 2011. By comparing the analysis results after MIP 
calibration, such as the visible energy, obtained with the MIP calibrations of the two independent 
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Table 4: Summary of uncertainties on the average visible energy and the average energy per layer 
due to the MIP calibration factors. 


Measuremenf 

Source 

Tofal 

Commenfs 

Average visible 
energy 

2% from correlafed sysfemafic 
uncerfainfies (til) 

±2% 


Average energy 

2% from correlafed sysfemafic 

±3.3% 

Added in quadrafure 

per layer 

uncerfainfies and 2.6% from 

uncorrelafed sfafisfical uncerfainfies 


due fo independenf 
origins 


periods, the systematic uncertainty of the MIP calibration can be extracted. As systematic uncer¬ 
tainty we use the width of the difference in the analysis results obtained with the two independent 
calibrations. This yields a variation of 2.6%, which is uncorrelated from one cell to the other*. 
This uncertainty has a negligible impact on the total energy sum, because the effect is reduced by a 
factor s/N, where N is the number of cells contributing to the measurement. However, its effect on 
the longitudinal shower profiles is larger, because the energy deposition in an individual layer can 
be dominated by few cells only, as in the case of electromagnetic showers. The uncertainties due 
to the MIP calibration factors are summarised in table 0. 

The MIP calibration factors are corrected for temperature variations as discussed in section 
using MIP temperature coefficients obtained by means of a linear fit. The average uncertainty on 
the MIP temperature coefficients is 2.3%. Modifying the MIP temperature correction within this 
uncertainty results in a negligible change in the average visible energy. 

Similar to the MIP calibration factors, the SiPM gain factors are also corrected for temperature 
variations. Uncertainties on the temperature correction of the SiPM gain are negligible, first, be¬ 
cause the gain variation itself is small — 1.7%/°C as compared to —4.5%/°C for the MIP value, and 
second, as they enter only via the saturation corrections, where uncertainties on the total number 
of pixels enter in a similar way and are dominant. 

6.3 Stability of the detector response in time 

The test beam experiments of the W-AHCAL at the CERN SPS were performed in two separate 
data taking periods in summer and autumn separated by five monfhs. We observed fhaf fhe defecfor 
response for dafa sefs af fhe same beam energy in fhese fwo periods was slighfly alfered, in par¬ 
ticular fhe dafa from negafive parficles were slighfly shifted with respect to the data from positive 
particles, independently of beam momentum. This cannot be explained by changes in the beam 
properties, consequently we attribute it to a detector instability. The origin and full time depen¬ 
dence are unknown, since corrections for temperature variation have been applied, and since the 
shift is not reproduced in the Monte Carlo simulations, which do account, for example, for changes 
in the beam position, beam profile and dead cells. The magnifude of fhe variation is beyond fhe 

*The uncertainty of 2.6% includes a contribution from the uncorrelated cell-to-cell variations of the MIP temperature 
coefficient, since we use a representative average coefficient per layer to correct for the temperature difference between 
the two periods. Hence, the overall uncertainty can he larger than the statistical uncertainty extracted from the MIP fit. 
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Figure 9: Histogram of the ratio of the visible energy per run over the visible energy for the full data 
set at one beam energy for (left panel) and (right panel). The relative number of events per 
run with respect to the full data sets at one beam momentum is used as weights for the histogram. 


corresponding uncertainties. 

Therefore we adopted the following strategy to quantify the stability of the detector response in 
time: we study the variation of the average visible energy deposited in the calorimeter measured in 
one run, (^vis, l■un)^ divided by the average visible energy measured using all runs at the same beam 
energy, (^vis, all at one energy ), including positive and negative polarity data, (^vis, run) is obtained in 
all cases with a Gaussian fit function performed in the central region of the energy sum peak de¬ 
fined by mean ±1.5 a. 

The obfained disfribufion of fhe W-AHCAL response including runs wifh posifive and negafive po- 
larify are shown in figure ^ for e* from 15 GeV fo 40 GeV (leff) and from 25 GeV fo 150 GeV 
(righf). The corresponding disfribufions in simulations exhibif only variafions of a few per mill. 

In fhe following, fhe sysfemafic uncerfainfy due fo fhe observed variafions of fhe calorimefer re¬ 
sponse is quanfified based on fhe sfandard deviafion of fhese disfribufion divided by fhe mean. 
For bofh e+ and hadrons, fhe resulfing sysfemafic uncerfainfy due fo fhe insfabilify of fhe defecfor 
response is 1.8%. The uncerfainlies are applied fo dafa. 

6.4 Shower start layer reconstruction 

As discussed in section P~2| , hadronic events are selected only if the shower starts within the first 
three W-AHCAL layers in order to reduce the effect of leakage. The observables used in this paper 
are tested for all beam momenta with data and simulation for two different shower start selections: 
shower starts in layers 1-3 and in layers 2-4. Differences of a few percent are observed. In general, 
the simulation describes these effects to better than 1 % except for the energy resolution where the 
difference is 2.1%. The observed disagreement between data and simulation is used as systematic 
uncertainty of the shower stait selection. The systematic uncertainty is applied to data. 
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6.5 Selection purity 


The analysed hadron data sets contain mixed muon, pion, kaon, and proton events. Events from 
a given particle type are selected using Cherenkov threshold counters and information based on 
shower shapes. While the muon contamination is efficiently reduced using calorimeter informa¬ 
tion, the distinction between pions, kaons, and protons is solely based on Cherenkov counters. The 
purity of the selected hadron samples was determined in [HBl . For this analysis, only runs which 
have a purity between 85% and 100% for a given hadron type are considered. An analysis of shower 
shapes in the Fe-AHCAL showed that for a particle purity above 85% the systematic uncer¬ 
tainty on the shower shapes is at most 1.5% on the shower radius and 1% on other observables. 
Since the profiles in iron and fungsfen are similar, one may assume fhaf fhe sysfemafic uncertain¬ 
ties due fo fhe parficle impurify are negligible in comparison fo fhe already discussed sysfemafic 
uncerfainfies. Therefore no sysfemafic uncerfainfies due fo possible confaminafions are considered. 

6.6 Systematic uncertainties of the simulations 

Optical cross-talk between tiles Due to the imperfect reflective coating of the scintillator tiles, 
light may leak between neighbouring calorimeter cells. This is taken into consideration in the 
simulation using the cross-talk factor, which is the fraction of light for which leakage will appear 
in neighbouring cells. Measurements of the cross-talk found values of 2.5% flT^ . and between 
3.3% and 4.6% lETH of the total energy per edge for the 3x3 cm^ cells. 

In the simulation, the cross-talk factor is relevant as it sets the energy scale. The higher this factor, 
the higher the final energy is on the MIP scale, since unlike for muon calibration events, in showers 
most of the signal leaks into neighbouring tiles which already have signals above threshold. 

To account for the imperfect knowledge of the cross-talk, simulations of positron- and hadron- 
induced showers are performed with cross-talk values corresponding to the lowest and the highest 
cross-talk measurements of 2.5% and 4.6%. The difference in the response divided by vT2 is used 
as an estimate for the systematic uncertainty due to the cross-talk. For both positrons and hadrons, 
a systematic uncertainty of 2.7% was found for the average visible energy. For the radial energy 
profiles studied in hadron-induced showers, a maximum systematic uncertainty due to cross-talk of 
4.2% was found. 

Timing cut As introduced in section ^ a timing cut of 150 ns is applied on simulated hits in order 
to emulate the signal shaping time of the CAFICE readout electronics. Due to differences in the 
SiPM production series used for the AHCAF, cell-by-cell differences in the signal shaping time 
of up to 15 ns were observed when optimising the readout time window of single cells. For the 
Fe-AHCAF, it was found that changing the timing cut within this observed variations of 15 ns has 
very small effects (at the per mill level) on the reconstructed hadron shower energy flTTIl . Because 
tungsten is used as absorber, a rather large fraction of the signal is produced by late spallation 
neutrons flia. While most of the energy is deposited within the first nanoseconds, there are also 
significant late energy deposits. To estimate the effects due to the timing cut, variations of ±30 ns 
were considered, resulting in variations of the measured energy in the simulation of ±0.3%, which 
is negligible. 
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Table 5: Systematic uncertainties on the W-AHCAL average visible energy {Eyis) for the beam 
momentum range of 15 GeV—40 GeV for and 25 GeV-150 GeV for K^, and protons. 


Source 


— Systematic uncertainty on (^vis) — 
for e+ (%) for 7i+, K+, and protons (%) 


Assigned to 


SiPM saturation scaling 1.4-3.0 

MIP constants 2.0 

Detector stability 1.8 

Shower start 

Cross-talk 2.7 


0.4-1.5 
2.0 
1.8 
0.1 
2.7 


data 

data 

data 

data 

MC 


Quadratic sum for data 3.1-4.0 

MC 2.7 


2.7-3.1 
2.7 


200 



400 600 800 

Energy sum (MIP) 



800 1000 1200 1400 1600 

Energy sum (MIP) 


Figure 10: Comparisons of the energy sum distribution of 15 GeV (left) and 40 GeV (right) data 

with simulation. 


6.7 Summary of uncertainties 

In this paper, an extensive study of systematic uncertainties of data and simulations is performed 
and summarised for contributions to the average visible energy (Syis) in table ^ The total uncer¬ 
tainty for a given data set is obtained by adding all relevant contributions in quadrature. 

7. Analysis of positron data 

In the following section, the analysis of positron-induced showers in the W-AHCAL is discussed. 
As Geant4 simulations of electromagnetic showers are well established, the analysis of the e+ 
showers permits a validation of the detector calibration and implementation of the detector simula¬ 
tion. Due to the limited penetration depth of positrons this is only possible for the first ~10 layers 
of the hadronic calorimeter. In addition, the W-AHCAL response to positrons is compared to the 
hadron responses, discussed in section ^ to study the degree of compensation of the calorimeter. 
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Figure 11: Comparisons of the longitudinal energy profiles of 15 GeV (left) and 40GeV (right) e^ 
data with simulation. Here and in the following figures, the uncertainty bands include statistical 
and systematic uncertainties while the data points only show statistical uncertainties. 


Figures |T^ and [H] show the energy sum distributions and the longitudinal profiles (i.e., the distribu¬ 
tion of the average energy deposited in a given calorimeter layer as a function of the layer number) 
for 15 GeV and 40 GeV positrons for data and simulations. The energy sum distributions of data 
and simulation agree reasonably well with each other. In the longitudinal profiles, discrepancies of 
up to 15% between data and simulation are observed. These are assumed to be due to missing in¬ 
dividual measurements of the SiPM saturation scaling factors, an imperfect description of material 
in the beam line, and uncertainties in Geant4 modelling. 

The average e+ visible energy (Fvis) and the width of the energy sum distribution Oe divided by the 
visible energy, i.e. the energy resolution, shown in figure |^, are obtained by fitting the e+ energy 
sum distributions to Gaussian functions. The fit has been performed in the central region of the 
energy sum peak defined by the mean ± 1.5 a. The average visible energy is given by the mean of 
the Gaussian fit function and the width of the energy sum distribution is given by the sigma of the 
Gaussian fit function. 

A minimum is used for fitting the calorimeter energy distribution and their energy dependence. 
The uncertainties of the parameters are obtained from the covariance matrix after the function 
has been normalised at the fit minimum to x^m = NDF (number of degrees of freedom) by scaling 

the uncertainties by the factor ^/zLnomuOised/NDF 

In the left panel of figure |T^ fhe average visible energy is shown as a function of the beam momen- 

^This method may incur problems if applied blindly Itfdll : here it is employed with due attention both for fitting the 
energy sum distribution and the energy resolution parameters a and b, listed in tables B0 and^ 
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Figure 12: Average visible energy (left) and energy resolution (right) as a function of the beam 
momentum for positron-induced showers for data and simulation. The black lines indicate fits 
to the data points without constraints. The bands show the overall uncertainty. For the energy 
resolution, re-analysed results from the corresponding PS data analysis i W^l are shown, as well as 
results from simulations with and without (default) detector instability (see text for definition). 


turn both for data and simulation. The systematic uncertainties are evaluated according to table ^ 
The experimental data for the visible energy agree well with the Monte Carlo simulations within 
the uncertainties, the deviations being less then 1%. The increase of (^vis) with beam momentum 
is in good agreement with a straight line fit shown for data as a solid line. The extrapolation of this 
fit to pBeam = OGcV is shown as a dashed line. 

The W-AHCAL e+ energy resolution is displayed in the right panel of figure |T^ To better constrain 
the energy resolution fit results, the low energy data from the PS test beam campaign, described 
in m, are included using the SPS selection cuts described in section M As discussed in section 
^ the simulation includes the effects of intrinsic calorimeter fluctuations as well as instrumental 
effects such as photo-electron statistics and noise. The detector instability estimated from the run- 
per-run variation of the calorimeter response in data (section |53| ), however, is not included in the 
simulation. This is referred to as default simulation in figure 0 In order to achieve a more re¬ 
alistic modelling of the data in terms of the energy resolution, the detector instability needs to be 
accounted for in the simulated energy resolution. This can be achieved by adding the observed 
effect of the detector instability on (^vis) from data in quadrature to the simulated positron energy 
resolution. After taking the detector instability into account, the energy resolution measured for 
positrons is reproduced by the simulations within uncertainties. 


- 19- 






















Table 6: Parameters of the energy resolution fits for data and simulation using beam momenta 
from 1 GeV to 40GeV The simulated results are obtained after including the detector instability 
measured in data. 


Parameter 

Data 

Simulation 

a (% • VGeV) 

29.5 ±0.4 

28.7 ±0.5 

b{%) 

L1±L1 

1.7 ±0.6 

c (GeV) 

0.035 

0.035 


The combined data from the PS and the SPS test beam campaigns are fitted to 


Oe a c 


(7.1) 


where a is the stochastic term and b is the constant term, which is dominated by the stability of 
the calibration; c is the noise term fixed to the noise in the fiducial volume of the e+ showers 
corresponding to 35 MeV, measured using dedicated random trigger events. The resulting energy 
resolution fit for data is displayed as a solid black line in figure |T^. The parameters of the positron 
energy resolution fits for data and simulation including the detector instability are listed in table ^ 
The indicated uncertainties include statistical as well as systematic uncertainties as introduced in 
section ^ The stochastic and constant terms derived from data and Monte Carlo agree within the 
uncertainties. The W-AHCAL energy resolution fit results are in agreement with those found using 
only the PS data sets IIT^ . The observed e+ energy resolution is slightly larger than that found 
for the Fe-AHCAL with a stochastic term of approximately 22% flT^ with beam momenta from 
10 GeV to 50 GeV. The larger value of a for the W-AHCAL in comparison to the Fe-AHCAL 
is expected due to the larger number of Xq per layer in tungsten (2.8 Xq) in comparison to steel 
(L2Xo). 

In summary, within uncertainties, the electromagnetic data agree well with the simulation for re¬ 
sponse, linearity, and energy resolution. For the longitudinal energy profiles layer-by-layer devia¬ 
tions of up to 15% are observed. 


8. Analysis of hadron data 

In the following section, the analysis of hadron-induced showers in the W-AHCAL is discussed 
and the results are compared to Geant4 simulations. First, the properties of pion-induced showers 
are studied, followed by shower properties for protons and kaons. 

For brevity, shower properties describing the evolution of hadron showers are presented only for 
pion-induced showers. The shower evolution of pion, kaon, and proton showers are consistent to 
within 20%. The agreement between data and simulation is at the same level for all three particle 
types. 

In order to take particles’ mass and decay properties into account, the observables are presented as 
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Figure 13: Energy sum distributions ofn^ events at 25 GeV (left) and 150GeV (right). The data 
are compared to selected Geantd physics lists. 


a function of the available hadron shower energy, £'available^ given by 

^'available = \/«^Lson+PLm (^-1) 

^available = "^Lryon + Pbcam “ "^baryon for prOtOnS, (8.2) 

where mmeson and mbaiyon are the masses of the corresponding meson or baryon and pbeam is the 
beam momentum. 


8.1 The pion data 

In the following, the calorimeter response, the energy resolution, and variables describing longitu¬ 
dinal and radial shower development of pion-induced showers are studied and compared to Monte 
Carlo simulations. The purity of the pion samples selected based on Cherenkov threshold counters 
is between 93% and 100% iHTIl . 


8.1.1 Calorimeter response and energy resolution 

In figure |T^, the energy sum distribution for 25 GeV and 150 GeV events are compared to 
the Geant4 physics lists QGSP_BERT_HP and FTFP_BERT_HP. Both physics lists describe the 
shapes fairly well. In the 150 GeV distributions of data and simulations, a tail towards lower 
energies indicates leakage effects in the calorimeter as discussed in section P~2| . 

The average visible energy (Fvis) deposited in the calorimeter by showers is estimated using 
Gaussian fits to the energy sum distributions as described in section ^ (Fvis) is shown in the range 
from 25 GeV to 150 GeV as a function of the available energy in figure |I^ (leff) where fhe available 
energy for 7t+ showers is given by equafion 8.1. The solid line indicates a linear fil fo fhe dafa and 
shows fhaf fhe response for pion-inifiafed showers is linear wifhin fhe uncerfainfies. The dafa are 
compared fo Geanl4 physics lisfs. A good agreemenf of ±2% befween dafa and QGSP_BERT_HP 
and FTFP_BERT_HP is observed which is covered by fhe sysfemafic uncerfainfies of fhe mea- 
suremenf. Af beam momenfa below 60 GeV, however, bofh Geanf4 physics lisfs fend fo slighfly 
underesfimafe fhe dafa. 
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Figure 14: Average visible energy as a function of the available energy (left) and energy resolution 
(right) for induced showers. For the energy resolution, data points from the corresponding 
PS data analysis /E2/ shown. The black lines indicate fits to the data points. The data are 
compared to selected Geantd physics lists. The bands show the overall uncertainties. 

Table 7: Parameters of the energy resolution fits for data and simulations using beam momenta 
from 3 GeV to ISOGeV. The simulated results are obtained after including the detector instability 
measured in data. 


Parameter 

Data 

QGSP_BERT_HP 

PTPP_BERT_HP 

a (% • VGeV) 

57.9±1.1 

51.1±2.8 

54.6 ±2.0 

b(%) 

4.6 ±0.4 

3.9±0.7 

3.8 ±0.5 

c (GeV) 

0.065 

0.065 

0.065 


The relative energy resolution for the 7t+ showers in the energy range from 3 GeV to ISOGeV is 
presented in the right panel of figure |I^ The simulation includes the detector instability, using the 
procedure established in section 0 and validated with positrons. The PS data were selected based 
on the new selection cuts described in section In the combination of the PS and SPS data for 
the energy resolution fit, the difference of the response between 30 and 38 layers has been taken 
into account by introducing a systematic uncertainty of 2.0% in the PS energy range. In the SPS 
beam momentum range, except at 150 GeV, the Monte Carlo simulations underestimate Oe/{E y{f) 
by 3%-12% for FTFP_BERT_HP and by 10%-15% for QGSP_BERT_HP. 

The energy resolutions for data and simulations are fitted using the energy resolution function 
defined in equation m The noise term c has been fixed to 65 MeV, corresponding to the electronic 
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noise in the full calorimeter volume. The resulting parameters of the energy resolution fits for 
all data sets are listed in table |7|. The stochastic term found for both physics lists underestimates 
the value found in data. The fit results of the stochastic term for the experimental data are slightly 
lower than those found in the corresponding PS data analysis WHU . This is due to the fact that for 
the PS study, the Ge / (£'vis) was estimated using the standard deviation and the mean of the energy 
sum distribution. These are more affected by the tails of the distribution than the parameters of 
the Gaussian fit functions which are used here for both PS and SPS data. The energy resolution fit 
results of the W-AHCAL agree with the corresponding fit for Fe-AHCAL combined with a steel 
tail catcher of in total 12 Aj obtained with beam momenta from 10 GeV to 80 GeV, with a stochastic 
term of approximately 58% \nM . 

8.1.2 Spatial development 

In order to assess the accuracy of the simulation of the spatial development of hadronic showers, we 
compare, for data and simulation results, the shower development along the z-axis (longitudinal) 
and in the xy-plane (radial). 

Longitudinal shower development. The development of hadron showers along the beam direc¬ 
tion (z-axis) is shown in longitudinal profiles, using as a reference fhe measured shower sfarf. Thus 
llucfualions of fhe shower sfarf are disenfangled from fhe infrinsic longifudinal shower develop- 
menf. The layer in which fhe shower sfarfed is identified using fhe clusfering algorifhm described 
in nTEH . We found fhaf dafa and Monfe Carlo af all beam momenfa agree in ferms of fhe recon- 
sfrucfed, exponentially falling shower sfait layer disfribufion wifhin 2% for shower sfarfs in fhe 
mosf relevanf firsf 15 layers. 

Longifudinal profiles mighf be particularly sensifive fo showers sfarting close fo fhe enfrance of fhe 
calorimefer. Therefore, an alternative even! selection is used here. Only evenfs sfarting in layer four 
or lafer are considered in order fo ensure fhaf fhe showers sfarf in fhe bulk of fhe calorimefer. Uncer- 
fainfies due fo layer-by-layer flucfuafions of fhe calorimefer response generafed by uncerfainfies in 
fhe calibrafion and safurafion effecfs are faken info accounf, in addifion fo fhe sysfemafic uncertain¬ 
ties discussed in section ^ The uncerfainfy due fo layer-by-layer flucluations is estimated by fhe 
difference befween fhe effecf observed in dafa and simulafions following fhe procedure described 

in UTEH . 

A comparison of fhe longifudinal profile in dafa and Monte Carlo for 25 GeV and 150 GeV 
evenfs is shown in figure [T^ Af bofh momenfa, QGSP_BERT_HP overesfimafes fhe energy depo- 
sifion in fhe firsf four layers which is compensated by an underesfimafe of fhe energy deposition in 
fhe subsequenf pari of fhe shower, while fhe opposife behaviour is observed for FTFP_BERT_HP 
af 25 GeV. FTFP_BERT_HP gives a good overall description of fhe longifudinal energy profile af 
150 GeV. 

A second characterislic parameter of fhe shower developmenf along fhe z-axis is fhe energy weighted 
cenfre-of-gravily in fhe z-direclion, Zcog^ as defined in equation P~T|. The Zcog disfribufion is shown in 
figure 1^ (leff) for 100 GeV evenfs. In fhe righl panel of figure fhe average cenfre-of-gravify 
in fhe z-direclion (zcog) > eslimaled as fhe mean of fhe Zcog disfribufion, is shown as a function of 
fhe available energy and compared fo fhe Geanf4 simulafions. FTFP_BERT_HP agrees wifh fhe 
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Figure 15: Longitudinal profiles from shower start for 25 GeV (left) and 150GeV (right) 7t+. 
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Figure 16: Centre-of-gravity distribution in the z-directionforpions of 100 GeV (left). Dependence 
of the average centre-of-gravity in the z-direction on the available energy (right). 


experimental data within 2% while QGSP_BERT_HP shows a 5% earlier centre-of-gravity in the 
z-direction than observed in data. 
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Figure 17: The radial profile of 25 GeV (left) and 150GeV (right). Here, the is visible 
energy {Eyifi normalised to the ring area. 


Radial shower development. The radial energy profile is given by the distribution of the energy 
density as a function of the radial distance to the shower’s z-axis. Here, the radial energy density, 
{E*ifi, is defined as the average energy sum in rings in the projection along the z-axis normalised 
to the ring’s surface. The radial distance to the shower centre is defined as 


n 



Y + iyi-ycogY, 


(8.3) 


where Xi (yi) is the x (y) position of the cell /, and Xcog CVcog) is the centre-of-gravity in x (y) 


•^COg — 


LjEj-Xj 

LiEi 


and jcog 


LiEi-yi 

HiEi 


(8.4) 


and Ei is the cell energy. 

Example radial profiles are displayed in figure for 7t+ at 25 GeV and 150 GeV. The bin width 
of the histogram’s horizontal axis corresponds to the dimension of the smallest W-AHCAL tile, 
i.e. 3x3 cm^. At both beam momenta, the Geant4 simulations predict a higher density in the core 
of the shower than observed in data. This effect decreases with increasing beam momentum. The 
better agreement with data is observed for FTFP_BERT_HP at 150 GeV. 
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Figure 18: Energy-weighted shower radius for pions of 100 GeV (left). Dependence of the average 
shower radius on the available energy (right). 


To further study the shower development in the radial direction, an energy-weighted shower radius 
is defined 


R = 


LjEj-rj 
LiEi ’ 


(8.5) 


where Et is the cell energy, and r,- is defined in equation |0| . A distribution of the energy-weighted 
shower radius is shown in the left panel of figure for 100 GeV tc+. The mean shower radius 
(/?), estimated as the mean of the R distribution, is shown as a function of the available energy 
in the right panel of figure |^. The better agreement with data is observed for FTFP_BERT_HP, 
within 7% or better, however both simulations significantly underestimate the shower radius at all 
energies. 
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Figure 19: Average visible energy (left) and energy resolution (right) for proton-induced (p) show¬ 
ers as a function of the available energy. For the energy resolution, data points from the corre¬ 
sponding PS data analysis /E2/ shown. The black lines indicate fits to the data points. The data 
are compared to selected Geantd physics lists. The bands show the overall uncertainties. 


8.2 The proton data 

The proton events are selected as described in Section The purity of the event samples selected 
based on Cherenkov threshold counters is between 85% and 100% UTOl . 

The average visible energy deposited in the calorimeter by proton-induced showers is shown as a 
function of available energy in the left panel of figure [T^ in the range from 25 GeV to 150 GeV. The 
data are compared to the selected Geant4 physics lists. 

The proton response in this energy range is linear within uncertainties. Data and both physics lists 
agree for all analysed energies within uncertainties. At beam momenta below 60 GeV, however, 
both Geant4 physics lists tend to slightly underestimate the data. 

The energy resolution for the proton data is presented in the right panel of figure |T^ Data in the PS 
beam momentum range from 4 GeV to 10 GeV nOl are included. As for the e+ and Tt+ data, the 
PS data were selected based on the new selection cuts described in section P~2| . The experimental 
energy resolution is up to 15% larger than predicted by the Monte Carlo simulations including the 
detector instability. FTFP_BERT_HP is closer to data (within 10%) than QGSP_BERT_HP, as 
already observed for the pion-induced showers. The parameters of the proton energy resolution fits 
using equation are listed in table |^. As also observed for the pion data, the Geant4 physics lists 
underestimate the stochastic term of the energy resolution fit. The fit results of the proton data are 
in agreement with the PS results when taking into account the different methods used to estimate 
Oe and (£'vis)^ as discussed in section |0| . 
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Table 8: Parameters of the proton energy resolution fits for data and simulations using beam mo¬ 
menta from 4 GeV to 150 GeV. The simulated results are obtained after including the detector in¬ 
stability measured in data. 


Paramefer 

Dafa 

QGSP_BERT_HP 

FTFP_BERT_HP 

a (% • VGeV) 

60.7 ±1.2 

56.6±1.7 

58.7 ±1.7 

b{%) 

4.3 ±0.4 

3.1±0.5 

3.7 ±0.5 

c (GeV) 

0.065 

0.065 

0.065 



LU 1000 2000 3000 

Energy sum (MIP) 

Figure 20: Energy sum distributions for 60 GeV K^. The data are compared to selected Geantd 
physics lists. 


8.3 The kaon data 

In the hadron beam, a small fraction of events are expected to have kaon-induced showers. Data 
corresponding to beam momenta of 50 GeV and 60 GeV were selected since the purity of the 
event samples selected based on Cherenkov threshold counters is 89% at 50 GeV and 99% at 
60 GeV llTHl . The energy sum distribution for 60 GeV K+ is compared to the simulation in fig¬ 
ure A good agreement with both Geant4 physics lists is observed. 

9. Comparison of the response for different particle types 

The average visible energy in e+, K+, and proton showers are compared in figure for dafa 

(leff) and Geanf4 simulations (righf). The upper parf of fhe figures shows fhe average visible en¬ 
ergy as a function of fhe available energy. The lower parf of fhe figures shows a comparison of 
fhe average visible energy fo fhe available energy given by ((^visinGev) — £'avaii)/£'avaii- For fhis, 
fhe average visible energy is converfed from MIP fo GeV using fhe e+ fif resulfs from 15 GeV fo 
40 GeV. The fif is exfrapolafed fo fhe full energy range of fhe dafa sefs. The bands show fhe overall 
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Figure 21: The response of the CALICE W-AHCAL to different particles as a function of the avail¬ 
able energy for data (left) and for QGSP_BERT_HP (right). The solid green line indicates a linear 
fit to the data. The dotted green line indicates the extrapolation of the line to 0 GeV and 150 GeV. 
The lower part of the figures shows the comparison of the average visible and the available energy 
given by {{EyisinCev) - Eavail) / Eavail, where {EvismCev) corresponds to (Eyis) converted from MIP 
to GeV based on the e^ fit parameters. The bands show the overall uncertainties. 


uncertainties. 

For data, the individual e+ points agree with the available energy within ±0.5%. Up to approxi¬ 
mately 60 GeV, the hadron results agree with the available energy within 2% which is well within 
the overall uncertainties of approximately 3.0%. For larger beam energies, the reconstructed en¬ 
ergies start to deviate from the available energy with values outside the systematic uncertainties 
of up to 3% and 4% at 150 GeV for protons and pions. Both longitudinal leakage and saturation 
effects contribute to the observed deviation from the linear behaviour at high beam energies. The 
calorimeter response to 7t+ for beam energies below 80 GeV is slightly higher than the response to 
protons, however, the data points agree within the uncertainties of the measurement. The responses 
of 71+ and protons approach each other again with increasing available energy. In conclusion, given 
the overall uncertainties in the measurements, the detector response is almost identical for all par¬ 
ticle types up to 60 GeV such that the W-AHCAL is found to be approximately compensating in 
this energy range. 

The observed behaviour is reproduced by the Geant4 simulations as shown in the right panel of 
figure A similar behaviour was also found for W-AHCAL test beam data and simulations of 
e+, 7I+, and proton beams in the energy range from 1 GeV to 10GeV HHIl . 
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10. Summary and conclusions 


We have presented analysis results of e+, 7r+, K+, and proton data with beam momenta from 
15 GeV to 150 GeV reeorded with the CALICE W-AHCAL prototype in 2011 at the CERN SPS. 
The measured ealorimeter response, energy resolution, and the longitudinal and radial shower 
development are eompared to predietions of the Geant4 physies lists QGSP_BERT_HP and 
PTPP_BERT_HP as implemented in Geant4 version 9.6.p02. The simulations show a reasonable 
overall agreement with experimental data. 

As eleetromagnetie showers are well understood, a eomparison of positron-indueed showers be¬ 
tween experimental data and simulations is used to validate the deteetor ealibration and the im¬ 
plementation of the deteetor simulation. Data and simulation agree with eaeh other within the 
uneertainties for the deteetor response and the energy resolution. Eor both e+ data and simulations, 
the deteetor response inereases linearly with the beam momentum. The W-AHCAE energy resolu¬ 
tion for e+ data and simulations results in a stoehastie term of approximately 29% and a eonstant 
term of approximately 1%. 

The longitudinal energy profiles for e+ data and simulations agree with eaeh other at the level of 
15%. This differenee eould be due to an uneertainty on the single eell ealibration faetors, imper- 
feet deseription of the material in the beam line in simulations, and to uneertainties in the Geant4 
modelling. 

When eomparing the positive hadron data to predietions of the Geant4 physies lists, similar eon- 
elusions ean be drawn for n^, K+, and protons. In the ease of the ealorimeter response to hadrons, 
a good agreement within ±2% between data and simulation is obtained for QGSP_BERT_HP and 
PTPP_BERT_HP. The W-AHCAE response is the same, within systematie uneertainties, for e+, 
7r+, K+, and protons for beam momenta up to 60 GeV above whieh leakage effeets in the hadron 
showers start to play an inereasingly important role. The level of agreement between the response 
for the partiele types studied is in line with the Geant4 predietions. 

Eor the hadronie energy resolution of the W-AHCAE, a stoehastie term of approximately 58% 
(61%) is found for tc+ (proton) data and a eonstant term of approximately 5% is found for both 7t+ 
and proton data. Eor both and proton events, the energy resolution has up to 5% lower values in 
ETEP_BERT_HP than observed in data, and up to 10% in QGSP_BERT_HP. 

The Geant4 physies lists reproduee the longitudinal shower development of hadron showers within 
15%. Eor the radial shower development, both physies lists prediet a higher density in the eore 
of the shower than observed in data, with a slightly better agreement, within 10%, observed for 
ETEP_BERT_HP. 

In eonelusion, we have shown that the Geant4 physies lists QGSP_BERT_HP and ETEP_BERT_HP 
reproduee average hadron shower properties reeorded with the highly granular CAEICE W-AHCAE 
at the pereent level, and spatial shower profiles af fhe 15% level or heller. The ETEP_BERT_HP 
physies lisl gives a slighfly heller agreemenl wilh dala for mosl observables and hadron lypes lhan 
QGSP_BERT_HP. The W-AHCAE response inereases linearly as a funelion of Ihe beam momen- 
lum, however al large beam momenla an onsel of leakage and saluralion effeels ean be observed, 
in bolh dala and simulations. The energy resolution oblained for Ihe W-AHCAE meels Ihe expee- 
lalions based on experienee from Ihe Ee-AHCAE. The e+ energy resolution of Ihe W-AHCAE has 
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a slightly larger stochastic term than the Fe-AHCAL energy resolution due to the larger number 
of Xq per layer in the W-AHCAL. The hadron energy resolutions for the steel and tungsten ab¬ 
sorbers agree with each other since the sampling in terms of nuclear interaction lengths is about 
the same for steel and tungsten. Based on the comparison between e+ and positive hadron data, the 
W-AHCAL is found to be compensating up to 60 GeV. 
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